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Abstract 

In this paper, we find all possible asymptotic behaviours of the solutions of 
the second Painleve equation y" = 2y 3 + xy + a as the parameter a —>■ oo in 
the local region x <C a 2//3 . We prove that these are asymptotic behaviours 
by finding explicit error bounds. Moreover, we show that they are connected 
and complete in the sense that they correspond to all possible values of initial 
data given at a point in the local region. 



1 Introduction 

The asymptotic analysis of the six nonlinear second-order ODEs (ordinary 
differential equations) known as the Painleve equations [|1| has been a subject 
of both classical @, [3| and modern [|], |5j interest. More recent interest was 
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revived by their importance in physical applications, particularly in statistical 
mechanics (see ||), and by their close relationship to soliton equations (see 

The solutions of the Painleve equations are, in general, highly transcen- 
dental functions known as the Painleve transcendents. An advantage of 
asymptotic description is that they give explicit behaviours of such tran- 
scendental solutions in terms of known classical functions. Such descriptions 
also give insight into the asymptotic behaviours of the solutions of soliton 



equations [§, [L0 



In this paper, we focus on the asymptotic analysis of the second Painleve 
equation (or Pjj) 

y" = 2y 3 + xy + a (1.1) 

in the limit as a — > oo. Pjj is well known to be a symmetry reduction of the 
modified Korteweg-deVries equation (and other soliton equations, see e.g. 
f|). In such a reduction, the parameter a arises as a constant of integration. 
Therefore, a large value of a is related to large values of initial data for the 
soliton equation. 

Although there have been many studies % g, [TO], [TT], § [TJ, § [TJ, |T§ of 
the Painleve equations in the limit as their independent variable approaches 
a singularity of the equation, the asymptotic limit as a parameter approaches 
infinity has not been widely studied, with two exceptions. Kawai, Takei, et 
al |ffBl [T7|, have studied the Painleve equations extensively in such lim- 



its through their associated linear equations, called isomonodromy problems 
TP| , Kitaev has studied the double asymptotic limit of the second 



Painleve equation as both x — > oo and a — > oo, again through an isomon- 
odromy problem. However, both studies appear to be restricted to formal 
special solutions. 

In this paper, we prove that in a region bounded by x <C a 2 / 3 the local 
asymptotic behaviour of the general solution as a — > oo is given by elliptic 
functions. We also show that there exist special solutions whose local asymp- 
totic behaviours are given by degenerate elliptic functions, i.e. singly periodic 
functions. These results are analogous to the results found by Boutroux 0, |3| 
for the Painleve equations in the limit as x — > oo. 

However, the domains of validity are quite different from those that are 
known for the limit x — > oo. In our case, these domains are full disks (punc- 
tured at poles and stationary points) of large radius, unlike the case x — > oo 
in which such domains are either bounded local patches or sectors of angular 
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width strictly less than 2tt. 

Our methods are direct and do not use the relation of Pjj to an isomon- 
odromy problem. Part of our motivation is to show that rigorous general 
asymptotic results can be obtained through natural direct methods without 
going through associated linear problems. 

However, we do rely on a characteristic property, i.e. the Painleve prop- 
erty, of the Painleve equations to obtain sharp error bounds. This is the 
property that all movable singularities of all solutions of each Painleve equa- 
tion are poles. In particular, we use the direct proof of the Painleve property 



of the Painleve equations given in p2|, to define local paths of integration 



and to estimate w on such paths. Error bounds can be obtained without 
assuming this property (see Boutroux|g). However, its usage allows us to 
obtain improved error estimates. 

The Painleve equations also possess many other deep properties. For 



example, their solutions can be written in terms of entire functions ||23|| , they 
possess exact special solutions in terms of classical special functions and 
transformations relating solutions to other solutions (see || for references). 
Asymptotics is a natural part of their study as nonlinear special functions. 

1.1 Asymptotic Regions 

The study of Pjj for a — > oo necessarily decomposes into three parts corre- 
sponding to three regions. First, consider a transformation of Pjj given by a 
dominant-balance analysis as a —>■ oo, namely 

y(x) = a l ^w{z) 

z = a l/3 x (1.2) 



which transforms Pjj to 



w" = 2w 3 + 1 + -w, (1.3) 

a 



where the primes denote z-derivatives. The three asymptotic regions are 
given by z <C a, z ~ a, and aCz. In this paper, we restrict ourselves to 
obtaining rigorous local asymptotic results in the first region. Note that this 
corresponds to x a 2//3 . 

These local results can be extended into the second region by using a 



nonlinear multiple-scales approach developed by Joshi and Kruskal [ I I 
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for the first and second Painleve equations. The results give insight into a 
well known asymptotic limit of Pjj. Consider |] the following transformation 

of I',, 

* - jo 
y(x) = eY(X) + 1 



e 5 



4 

TTF 



e 1 

In the limit as e — > 0, this leads to the first Painleve equation (Pj) 

Y" = QY 2 + X. (1.4) 

It will be shown in a subsequent paper that this limit occurs at isolated 
points in the second region. 

The method of Joshi and Kruskal was developed for the limit x —>■ oo. 
The usual asymptotic studies of the Painleve equations in the limit x — > oo 
apply here in the third region given by ft 2 / 3 <C x. To see this, consider the 
Boutroux transformation of Pn that makes the maximal dominant balance 
as x — > oo explicit, namely 

y(x) = \fxu(t) 

t = -x 3 / 2 (1.5) 

3 v ' 

under which Pjj becomes 

u u = 2u z + u + —^ + — , (1.6) 

where a = 2a/3. (See for a deduction of Boutroux transformations based 
on maximal dominant balances.) The asymptotic analysis of this equation 
is usually carried out for t — > oo. It remains valid so long as a/t 1, i.e. 

ft 2/3 < x. 

1.2 Main Results 

In section 2, we study the generic solutions of Pjj. These are solutions defined 
by a large open set of initial data given at a point Zq in the local region. Such 
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solutions are shown to have asymptotic behaviours given by elliptic functions 
as a — > oo. In section 3, we study limits of the generic solutions and their 
asymptotic behaviours. The union of the results of Sections 2 and 3 cover 
all solutions defined by all possible initial data at zq. 

Moreover, the domain of validity of the generic locally valid elliptic- 
function-behaviour consists of the entire local domain z <C a (punctured 
at poles and stationary points of the solution). Similarly, its degenerate 
limits are also valid in large subdomains of this region. 

The solution in the nearly degenerate case is described by an almost- 
singly-periodic behaviour, exhibiting widely spaced lines of poles. There are 
further limits which exhibit only one line of poles or no such line within the 
domain of validity. We also consider the large energy limit and show that 
the resultant behaviours are still given by (scaled) elliptic functions in large 
subdomains. 

Our results are given for both the solution and its first-derivative. By 
using Pjj, these may be extended to yield the asymptotic behaviours of any 
higher-order derivative of w(z). 

In summary, we find all possible asymptotic behaviours of true solutions 
in the limit as a — > oo in large domains lying in the region z <C a. We prove 
that these are asymptotic behaviours by finding explicit error bounds. We 
show that these behaviours are connected through limiting values taken by 
an energy-like parameter and are complete in the sense that they correspond 
to all possible values of initial data given at a point Zq with z <C a. 



1.3 Notation 

It should be noted that we use two symbols for an asymptotic balance be- 
tween two functions f(s) and g(s) as s — » oo (along some path in the complex 
s-plane). First, if the limit of f/g is unity then we write / « g. Second, if 
a nonzero limit exists but differs possibly from unity, we write / ~ g. Our 
remaining usage of asymptotic notation is standard and definitions may be 
found, for example, in |E3 . 



2 Generic Local Asymptotic Estimates 

In this section, we prove that in the region where z <C a, the general solution 
of Pn is asymptotic to Jacobi elliptic functions. In particular, we estimate 
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the size of the asymptotic correction term. Moreover, we show that the 
spacing between zeroes of the general solution of Pjj is closely approximated 
by a period of the locally valid elliptic function. It should be noted that we 
obtain an asymptotic result for both w'(z) and w(z). Note that we assume 
| at | > 1 throughout. 

First, we describe the local domain for our asymptotic analysis. 

Definition 2.1 Suppose e 7^ 0, z$ are given complex numbers and B is a 
positive real number. They are called appropriate if 



< 1, 



In el 



a 



< e 



< 



, B 
|e|,— < 2 e 

a 



(6 + 271-5)2^5 < 
For an appropriate set {e, Zq, B} 

Z := {z 
is called the local z-domain. 



- 5 / 7 llne| 



-2 



zol < B}, 



Remarks: 

1. The last inequality above is imposed to ensure that the error estimates 
obtained in Theorem 2.1 below is small. See further remarks following that 
theorem. 

2. Note that one choice of appropriate {e, z ,B} is given by e = l/y/a with 



a 



x l\ B 



a 



3. The above remark shows that zq and B do not have to be small. In par- 
ticular, the local domain may be quite large (but not unboundedly large) 
even if z is finite. 



Standard theorems [p^j show that for given bounded initial data w(z ), 
w'(z ) at a finite point z there exists a unique locally analytic solution w(z) 
of Pjj. It is well known that w(z) can be meromorphically extended ]22] to 
any bounded domain. In particular, w(z) is meromorphic in the local domain 
Z. 
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We will work with an integrated form of Pjj. Multiply Eqn( |1.3| ) by w' 
and integrate as though the small term zw/a were not there. Then, we get 

w' 2 = w 4 + 2w + 2E + 2 T — dz, (2.7) 

Jzo Oi 

where E is a constant of integration. To leading-order, this equation is solved 
by Jacobian elliptic functions. 



2.1 Elliptic Functions 

To see that elliptic functions arise, and to fix the notation, let 

P( w ) ■= w 4 + 2w + 2E. 

In general, P(w) has four distinct roots, which we will denote by Wi, i = 
1, . . . , 4, in the complex plane. The leading order equation 

w' 2 = w 4 + 2w + 2E, (2.8) 

can be integrated again to give 

f w dv , s 

\ ^= = z-z , (2.9) 

where rj = w(zq). The path of integration here lies on the Riemann sur- 



face defined by JP(v). (This consists of two copies of the complex plane 
punctured at the roots of P(v) and glued along lines connecting two pairs of 
roots.) 

There exist two linearly independent closed contours Cj, j — 1, 2, on this 
surface each enclosing two roots of P(v) in the f-plane. Paths connecting r\ 
to the same end-point w are deformable to each other modulo Cj. That is, 
deformations lead to Eqn (|2.9| ) where w remains unchanged while the right 
side changes by integer multiples of 

u j = I dV . (2.10) 



Further properties of the elliptic integrals uij, including differential equations 
they satisfy as functions of E, are given in Appendix A. 
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This shows that the solution w(z) (the inverse of the function z(w) defined 
by Eqn (|2.9| )) has two periods given by Uj. Moreover, it is straightforward to 
show that w(z) is meromorphic. That is, it is an elliptic function. 

One of the two periods of w becomes infinite (while the other remains 
finite) at special values of E. These are the values for which two roots of 
P(w) coincide. To find these, consider 

p'( w ) = Aw 3 + 2 = => w = d k := (-1/2) 1/3 , k = 1,2, 3. (2.11) 

If such points are also zeroes of P(w) then we must have 

o / 07 \ 1/3 

E = D k := --(-l/2f 3 = j , k = 1, 2, 3. (2.12) 

For each such value of E, P(w) has a double root. Note, however, that P{w) 
cannot have a triple root, because P"(w) = 12w 2 cannot vanish simultane- 
ously with P'(w). 

Suppose E = Dk, for some k = 1,2, or 3. Let p and a be the roots of 
P(w) distinct from dk- Then rewriting 

P(w) = (w - d k ) 2 (w - p)(w - a) 

in expanded form shows that 



Solving these equations, we get 

,,„ = -L{l ±i V2} 



Note that these cannot equal d k because P{w) cannot have a triple root. 
That is, there exists c > 0, such that \p — d k \ > c, \a — d k \ > c. 

Now suppose E = D k + S 2 , for 5 <C 1. Then setting P(w) = and using 



the fact that d k — p, d k — a are bounded below, it is straightforward (see |26 



for detailed proofs of similar results) to show that the two nearly coalescing 
roots Wj of P(w) satisfy 

wj = d k + 0(5), 

and that the other two roots are a finite distance away from d k . 
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Note also that if E is large, say \E\ = log 8, then the periods become 
small. To see this, consider the roots of P(w). The roots now have modulus 
at least equal to OdE] 1 ^ 4 ) and so 

u j = 0(\E\-* / %j = l,2- 

(Such asymptotic results for the complete elliptic integrals can be obtained 
straightforwardly from the differential equations derived in Appendix A.) 
Hence the poles of the elliptic function crowd together and the function itself 
becomes large even at points of holomorphy. 



2.2 Generic Solutions 

Our main result holds for solutions that are described by a large open and 
dense set of initial values that we call generic. 

Definition 2.2 Suppose e is an appropriate given number. Parameters n, n' 
are called generic if for 

2E = r,' 2 -r, 4 - 2r], 

they satisfy 

\e\ 1/7 < \n'\ < | In e| 2 , \n\ < | lne| 
2|e| 2/7 < \E-D k \ 

These supply the initial values for the solutions and ensure that their formal 
elliptic-function-type behaviours are not too close to degenerate ones. Initial 
values that violate this definition are studied in section 3. 

Definition 2.3 Suppose e, z Q are appropriate numbers. Let w(z) be a so- 
lution of Pjj defined by generic initial values n, n' at z , i.e. w(z ) = rj, 
w'(zo) = rj' . Then w(z) is called a generic solution. 

To get valid asymptotic results, we need to confine such a solution w(z) to 
domains that avoid its poles and stationary points. 

Definition 2.4 Given an appropriate e, define the w-domain by 
W := {w I \w\ < | lne|, \P(w)\ > |e| 1/7 } 

and z-domain by 

Z:={zeZ\ w{z) e W}. 
e is called small enough if Z is connected. 
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The direct proof [^2] of the Painleve property of Pjj shows that there exists 
an open disk of values e for which Z is connected. (Explicit lower bounds 
for A := l/|e|, may be found on p. 197 of P2]| .) Moreover, any two points in 
Z are connected by paths 7, of bounded length 2nB, such that their images 
T under w lie in W with length < 2ir\ lne|. 
Our major result is 

Theorem 2.1 Given an appropriate set {e, zq,B} where e is small enough, 
and generic initial values r], r{ , the generic solution w(z) satisfies 



w 



rl 



-P(w) 



< 



K , 



In el 



dv 



PM 



{z - z ) 



< 2vr/t J B|e| 6/7 |lne| 2 



for all z G Z, where n := Q + 2ttB. Moreover, if zq and z\o are two successive 
points in Z where w = rj, then for j = 1 or 2, 

\{z w - z ) - Uj\ < 27fKB\e\ 6/7 \ lne| 2 . 



Proof: Equation( [2.7D can be written 

(w'(z)f = P{w{z)) + 2S 

where 

S :— / —ww'dz 

J-y a 

where 7 is a path connecting z to z in Z. The assumptions give 



(2.13) 



B 

\a\ 



< 3|e|. 



Integration by parts shows that 

S 

Hence we get 



zw 2 


(z,w) 


"/ 


2a 


(20, v) 


2a J 



2\S\ < 6|e| I lne| + |e| • 2nB\ lne| 
= /tlel I lne| 2 
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as desired. Now integrate Eqn ([2.7|) once more by taking the square root and 
dividing by JP(w). We get 



w dv 



P(v) 



, L 25 
dz { H -- 

7 I P(W) 



1/2 



(2.14) 



Note that we fix a choice of branch of the square root here, so that as e — > 



we have the positive real branch. From Eqn (|2.14| ) we get 



u dv 



PM 



{z - zo) = L dz \[ 1+ p7w) 



1/2 



Note that 



25 



P(w) 



<de| 6 / 7 |lne| 2 . 



(2-15) 



(2-16) 



Let 



Then we have 



Q:= 1 + 



25 



1/2 



+ 1. 



Q 2 = 2Q 



25 
P(w) 



Q = 2 



2S 



QP(w)' 

The choice of branch ensures that Q is lower bounded by \Q\ > 1. Therefore, 
it follows that 



\Q-2\< 



2S 



P{u) 



1 



2S 
T(u) 



1 



< K |e| 6 / 7 |lne| 2 . 



Integration then gives the desired result. The separation z 10 — z is based 
on the same result with the path of integration T being a closed contour 
enclosing two roots of P{u). □ 



Remarks: 

1. The crucial element of this proof lies in the lower bound for P(w). The 
microscopic details of this lower bound can be changed without affecting 
the macroscopic form of the asymptotic result. For example, suppose the 
conditions on i], r]' are changed to allow E — Dk = S 2 where \S\ ~ |e| p,/7 , 
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p < 7. Then the same result as in the theorem holds for u>' 2 , and the second 
(and third) results of the theorem hold with the order of the error estimate 
changed to 

H^llnel 2 . 

2. It should be noted, however, that if 5 is small enough (but nonzero), then 
the incomplete elliptic integral on the left of the first result is more accu- 
rately written as the inverse of a hyperbolic function. See the next section 
for details. 

3. The third result of the theorem holds for any two successive points (in a 
direction given by ujj for a chosen j) at which w takes the same value. By 
taking e to be smaller, such points can be taken arbitrarily close to poles 
of w. Hence the result also holds for the spacing between two poles. 

4. Notice that the incomplete elliptic integral 

r w dv 
h J¥(v) 

is still convergent if the generic initial value r\ happens to be arbitrarily 
close to a zero of P(w). However, if it is so close to such a zero that the 
path intrudes into a "hole" of Z, then we need a slightly modified argument 
for the proof. In this case, evaluating Pjj at zq gives a finite nonzero lower- 
bound for |u/'(zo)|- Therefore, there exists an ascent path for \w'\ starting 
at Zq. Following this path (over a bounded distance) takes the solution to 
a point z\ where P(w) is lower-bounded by \e\ 1 ^ 7 ■ We can now start again 
at such a point and apply the theorem. Such arguments are used in more 
detail in section 3. 

5. In the proof above, B is not required to be finite, but it should be appropri- 
ate. As remarked after Definition |27I|, B may be quite large. For example, 
the error estimates of the theorem are still small for B m e -1 / 4 « a 1 / 8 . 

The last remark above shows that for any finite Zq and large but appropriate 
B, a generic solution w has the same elliptic-function-behaviour throughout 
the "local" domain Z. The error, however, does increase with B. If B 
increases beyond appropriate values, z encroaches on the second asymptotic 
region (i.e. z ~ a) which will form the subject of a subsequent study. 



12 



3 Limits of Generic Solutions 



In this section, we consider the possible asymptotic behaviours of generic 
solutions w(z) when the parameter E or initial values 77, 7/ become close to 
degenerate values or become large. 

Degenerate limits occur when E — Dk =: 5 2 , with \S\ <C 1, for some k — 1, 
2 or 3. We separate this into the following subcases. 

1. For some integer p > 1, |e| p / 7 < \5\ < (el 1 / 7 

2. |<5| < |e| p / 7 , for any integer p > 1 

3. For some = 1, 2 or 3, |t/| < le) 1 / 7 , \rj — < | e | 1 / 14 

Note that these cases overlap. The first subcase is studied in Subsection 3.1 
and the remaining subcases are studied in Subsection 3.2 below. 



The results of the direct proof of the Painleve property of Pjj |22| show 
that large initial values with finite E imply that the initial point z is close 
to a pole of the generic solution w(z). The proof is a trivial rewriting of the 



proof given in section 3 of |22| and so we do not reproduce it here. In this 
case, Theorem [2.1| shows that the pole of the Painleve transcendent nearest 
to zo lies close to a pole of the elliptic-function-type asymptotic behaviour of 
w(z). 

Therefore, the limit of interest, that violates the genericity condition, is 
when E becomes unbounded. (Equivalently, we may take rj, rj' large with 
7/ 7^ i] 4 .) In Subsection 3.3, we consider the large energy limit when \E\ > 
I In el 4 . 



3.1 Nearly Degenerate Solution 

The results of Theorem |2.1| continue to apply even if 5 is smaller than e 1 / 7 . 
(See Remark 1 following Theorem |2.1| .) However, for small enough 5, the 
elliptic- function- type behaviours of w(z) are more accurately represented by 
singly periodic behaviours. We show this here. 

From the results of Section 2.1, we have that if E = Dk + S 2 , where 8 <C 1, 

P(w) = (w- d k f(w -p)(w-a) + 25 2 . 

Recall dk = (— 1/2) 1 / 3 , and p, a equals one of 



24 
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Definition 3.1 Let {e,z ,B} be an appropriate set, k = 1, 2, or 3, and 

1 < p be a given number. Parameters rj, rf are called nearly degenerate if for 



2E = rjf 



12 



V 



2rj, E — D k =: 5 2 , 



they satisfy 



|V7 



< \r)'\ < | In ej 



1 1/14 



< \v — dk\ < | lne| 



|e| p/T < \S\ < \e\ 1/7 



Note that the condition on E implies r]', rj lie in a one-complex dimensional 
deleted neighbourhood which excludes the point rj' = 0, rj = d k , for a given 
k. This deleted neighbourhood of initial conditions gives rise to a solution 
that is closer to a singly periodic limit of the previous elliptic function. 

Definition 3.2 Suppose e, zq are appropriate numbers and rj, rj' are nearly 
degenerate parameters. Then the solution w(z) of Pjj given by w(zq) = rj, 
w'(z ) = rj' is called a nearly degenerate solution. 

Now the appropriate polynomial to consider is 

R(w) := (w - d k ) 2 {w - p)(w - a), 

and the W-domain becomes 

V = {w J \w\ < \ \ne\,\R(w)\ > e 1/7 }. 

Correspondingly, e is called small enough if the preimage Z$ of V in Z is 
connected. To be specific, we take k — 1, with real 

d 1 = -1/2 1 / 3 . 

The main result in this case is 

Theorem 3.1 Given an appropriate set {e,z ,B}, with n\e\ 5 ^ 7 \ In e| 2 < 1, 
where e is small enough, and nearly degenerate parameters r\, rj , the nearly 
degenerate solution satisfies 

• /2 R(w) <3|e| 2 / 7 , 



w ' 



dv 



v - d^Mv- p){v - a) 



~(z- z ) 



< 67iB\e\ l/7 , 



for all z G Z$. 
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Proof: The proof is almost identical to that of Theorem [2.1| with a modified 
small term. Note that 

w' 2 = R(w) + 25 2 + 2S, 
where S is as defined in Eqn Q2.13|) . Therefore, since we still have 

2\S\ < /t|e||lne| 2 

whereas 

\5\ 2 < |e| 2 / 7 , 

we get the first result. Integrating once more, we get 



/ 



w dv 



'Rip) 

By assumption, we have 
and 

So we have 



, T 28 2 2S 

dz<l + -^7— r + 



1/2 



25 2 



R(w) 



R(w) R(w) 
< 2H 1 / 7 



(3.17) 



2S 



R(w) 
25 2 



< K |e| 6 / 7 |lne| 2 . 

■f^r <3|r|' 7 . 



R(w) R(w) 

The remainder of the proof of Theorem [2.1| can be repeated with {25 2 + 
2S} / R{w) in place of 2S/P(w) and R(w) in place of P(w) to get the desired 
result. □ 

Note that we now require a smaller B to get asymptotic validity. Namely, 
we need QnBe 1 ^ 7 « 1. For example, 6ttB < |e|~ 1/l4 would be sufficient. 

The interesting difference from the generic case arises because the integral 
on the left can be evaluated explicitly in terms of hyperbolic functions. To 
see this, make the change of variables: 



/ 



(v-Pl f f + f- 1 
2 V 2 

V3g-V2i 



(a + p) 



(3.18) 
(3.19) 
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Then the integral becomes 



dv 2 . f dg 



v (v-d^Mv- p)(v-<j) c?iV6 J g 2 + i 



Hence we get 



g = itanh {^±{z - zq) + a + 0(e^ 7 )) , (3.20) 



where a = tanh _1 (g(z )). Note that this behaviour is singly periodic and 
possess a line (rather than a lattice) of poles. In a direction orthogonal to 
the line, the behaviour becomes asymptotic to constants. 

There is an alternative change of variables. Consider the transformation: 

v = d l + u (3.21) 
u = -1/F (3.22) 

' G+^- (3.23) 



'18d? 3rf i 

Then the integral becomes 

rw-h du 1 r dG 



'v-di u^u 2 + Ad lU + 6dj VQdJ VG 2 + 1 ' 
In this case, we get 

G = sinh [VQd^z - zq) + b + 0(e 1/7 )} , (3.24) 

where sinh6 = G(z ) = —\/2 + \^18di/(d 1 — rj). The behaviours indicated 
by Eqns (pr20| ) and ( |3.24j ) are clearly related and represent the same nearly 
degenerate solution in the domain Z$- Note that one line of poles is given 
approximately by 

-b + 2nni 

20 + V6rf, <3 ' 25) 

for integer n. 

In Appendix A, we show that one of the periods uj of the generic elliptic- 
function behaviour becomes logarithmically large and one remains finite as 
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E — > for some = 1, 2, 3. For the nearly degenerate case this shows that 
for some j 

ujj ~ In | E — Dk | ~ | In e | . 

For sufficiently small e, this is clearly smaller than B (e.g. B ~ |e| -1 ^ 14 ) for 
which the asymptotic results of Theorem |3.1| are valid. Therefore, the results 
above suggest that more than one line of poles may be visible in the domain 
Z$. Between the lines of poles, the local behaviour is still described by Eqn 

(EH )- 



3.2 Degenerate Solution 

Here we consider the case when 5 is arbitrarily small or when initial values 
are arbitrarily close to degenerate values. 

Definition 3.3 Let {e, zq, B} be an appropriate set, and, for some k = 
1, 2, 3, let i], rj be given parameters. 

• // they satisfy 

2E : = r/ 2 - if - 2rj = 2D k + 25 2 , 
where \5\ < |e| p / 7 for all integers p > 1, and 

\e\ 1/7 < \r}'\ < | lne| 2 , \e\ l ' lA < \ V - d k \ < \ lne| 

then they are called almost degenerate parameters. 

• If, on the other hand, they satisfy 

|VI<H 1/7 , |r/-4|<|er/ 14 

then they are called degenerate parameters. 

To be concrete, we take k = 1. Assume that e, V, Z$ are as defined in the 
previous section, with a slight modification that \R\ > lej 1 / 7 /^ on Z$. 

Definition 3.4 Let e be small enough such that the domain Z$ is connected 
and let 7] v( be almost degenerate parameters. Then the solution w(z) of Pjj 
given by w(z ) = r\, w'(zq) = r/' is called an almost degenerate solution. 

The solution defined by almost degenerate parameters has asymptotic be- 
haviour that is given locally by a singly periodic function. In fact, the proof 
of the last subsection can be applied to show that the following theorem 
holds. 
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Theorem 3.2 Given an appropriate set {e, z ,B}, with /t|e| 5 / 7 | lne| 2 < \, 
where e is small enough, and almost degenerate parameters rj, rj' , the almost 
degenerate solution satisfies 



w' 2 - R(w) 



K | 

2 1 1 



In el 



dv 



>V (v - cZi) y/(v-p)(v-(7) 
for all z G Zg. 



< 27Tfi: J B|e| 6/7 |lne| 2 . 



Note that the domain of validity is now larger than in Theorem |0 

lei 5 / 7 . 



i.e. 



B 



may be as large as 1/ 1 

The almost degenerate case clearly also has representations given by 
Eqns( j3~20| ) and (|3.24|) in the domain Zg. However, now there may be at 
most one line of poles. The reason lies in the fact that the almost degenerate 
limit of an elliptic- function behaviour has a period that is possibly ^> | lne|. 
For example, suppose 5 ~ exp(— l/|e|). Then 



\n(E - D k 



Hence the distance between two lines of poles (or any two lines formed by 
interpolation through integer multiples of the finite period) may be much 
larger than appropriate values of B. 

Nevertheless, there must be at least one line of poles in the local domain, 
because the lower bound on rj — d\ implies |&o| < I hie|. Hence the sequence 
given by Eqn( p.25|) lies in Z§. 

The remaining case, apparently different to this, is the degenerate case. 
However, we can show that the almost degenerate solution becomes a degen- 
erate solution somewhere in its domain of validity. To show this, we travel on 
a descent path. Recall that a complex-valued function f(z) of one complex 
variable z possesses a descent path (for its modulus) starting a point Zq if 
f'(zo) 7^ 0. Along the steepest descent path, we have \df\ = —d\f\. 

Theorem 3.3 Let {e, zo, B} be an appropriate set with jej 1 / 4 < 2~ 1//3 . Let 
rj, rj' be almost degenerate parameters and let w be an almost degenerate 
solution. Then there exists a path T joining z to a point z 1 where \w — c4| 2 
and \w'\ become less than \e\ 1 / 7 , for some k = 1,2,3. 
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Proof: Case 1: Consider the case when \r}\ > |(4|/4. By assumption of 
appropriateness and upperbounds on rj, rj', we have 



V + z i]' 



a 



< 2|e||lne| 2 . 



Therefore, we can lower-bound 



w'" = 6w w' H — (w + zw') 



a 



(3.26) 



at (and by continuity near) z by 



\w 



\z )\ > |4| 2 | e | 1/7 /4. 



In other words, there exists a descent path for w" starting at z . We take T 
to be this descent path. 
Recall that in Pn 



and its first integral 
we have 



z 

—w 

a 



w " = 2w i + l + -w 

a 



w u = R(w) + 2S, 



<3|e||lne|, \S\ < «/2|e| | lne| 



(3.27) 
(3.28) 
(3.29) 



throughout Z$, and in particular on the part of the path T lying in Zg. 

Since w" is decreasing, this shows that 2w 3 + 1 must decrease along T, 
until it reaches a first point Z\ where \w — d k \ becomes smaller than (e] 1 ^ 14 , 
or where we reach the boundary of Zg. Similarly R(w) also decreases along 

r. 

At that point, the definitions of c4, p, and a show that 



R(w(zi)) 



1 1/7 



> 



16 ' 



while it remains larger than this bound earlier on T. Moreover, we can 
estimate its length. For conciseness, we rename v := w". By definition we 
get 

rzi 

z\ — z = / dz 

Jz 



r dv 
Jr v' 
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But from Eqns (13~2^ ), ( ggg ), and ( ggg ), we get 



u'| > 4|tu| VI > |4| 2 |er /14 /64 



if 



| e | 1/14 >32vr 2 rf 2 | e ||lne| 
where d is the length of V. We then get 



(3.30) 



d 



< 



< 



dv 

v v' 
\dv 



r \V 
64 



|4| 

64|er 1 / 14 
I4I 2 



d|«|/|e 



1/14 



(|u(2b)| - |«(«l)| 



This clearly satisfies Eqn (|3.30| ). Moreover, this shows that T lies within Z$. 
Case 2: Now consider the case when \rj\ < |rffc/4|. We have from 



that 



w " = 2w 6 + 1 + -w 

a 



\w"(z )\ > 1/2. 



Hence there exists a descent path for w' starting at zq. We take T to be 
this path. The decrease in \w'\ implies a corresponding decrease in 
Therefore, w must get closer to a root of R(w). Since all such roots have 
modulus > 1/2 1//3 , \w\ must increase. There exists a first point Zq on T where 
w has modulus equal to |<4/4|. That is, 



\w(zq) - 41 > 



3|4 



At such a point, we have 



\R(w)\ > 18|4 



by using the definitions of 4, P, and a. Using arguments similar to above 
we get 



I — Zq\ 



--:d< 



1 

2VU 
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We now use the results of Case 1 starting at z to extend r to z\. □ 

It is instructive to consider the limits of Theorems |3.1| , |3.2| , explicitly 
when 

^ \ w - 41 < H 1/14 , 0^\r]-d k \< (el 1 / 14 . 



We have 



w dv 



R{w) 

dv 



VQd k J v (w - 4)^1 + 2(w - 4)/(34) + (w- d k ) 2 /(Qd 2 k ) 

= ^J^^{ 1 - {w - iM+0((w - ik)2) } 

= -^-|ln((«; - d k )/( v - 4)) -(w- r/)/(34) + 0((w - 4) 2 )| 

Taking the inverse of the second result of Theorem |3.2| implies that 
w-d k = (r/-4)exp(V64(^-^o)) (l-(^-4)/(34) + 0(|e| 1/7 )) 
+ ^-(r/-4) 2 exp(2V64(^-^o))(l + 0(|e| 1 / 14 )). (3.31) 

It is clear that the exponential terms here will grow in a half-plane of direc- 
tions in Z$. For example, consider 4 = — 1/2 1 ' 3 . Then 

$l(z - zq) < 3?(4(z - zb)) > 0. 

In this half-plane, the solution grows away from the initial assumed value. 
By following an ascent path, it is straightforward to show that this behaviour 
matches with that given by Eqn( |3.24| ) in this half-plane. 

Conversely, if the solution has degenerate initial values, we can integrate 
Pjj to show that w is given by Eqn (|3.31|) in a local domain around zq. To 
prove this, we consider a domain given by 

U = {w | |e| 1/7 <\w\<\ lne|, \w - d k \ < e 1/14 }. 

Correspondingly, define e to be small enough if the preimage Z ofUinZ is 
connected. 



21 



Definition 3.5 Lett be small enough and letrj, 77' be degenerate parameters. 
Then the solution w(z) of Pjj given by w(z ) = n, w'(z Q ) = n' is called a 
degenerate solution. 

Then we have the following theorem. 

Theorem 3.4 Assume we are given an appropriate set {e, zo, B}, where e is 
small enough such that 

/t|e| 5/7 |lne| 2 < 1, |e| 1/7 < \d k \, 



and 



6|4| 2 + 4|4||e 



1/14 



+ e 



1/7 



<7\dl 



Furthermore, assume that degenerate parameters n 7^ dk, rj' 7^ 0, are given. 
Then there exists a number fi > such that the degenerate solution satisfies 



pi 

Jn 



w' 2 — 6dl(w — d k ) 2 
dv 

Qd 2 k (v - d k f + 25 2 



1/14 



for all z £ Zo- 



Proof: Again, the technique of the proof is almost the same as Theorem 2.1 



and we only highlight the important differences. After one integration, we 
get 



w 



R{w) + 25 2 + 2S, 



where by the given assumptions, we have 

2\5\ 2 = 2\E-D k \ < 8|4| 2 |e| 1/7 . 

S here is as before, with the same upperbound: 1251 < fc|e|| lne| 2 . To get the 
first result, we use the fact that 

R(w) = 6d 2 k (w - d k ) 2 + 4d k (w - d k f + (w - d k ) 4 , 

and the bounds oiimj-4 

Hence it is slightly more transparent to write 

w' 2 = 6d 2 k (w - d k ) 2 + 440 - d k f + (w- d k f + 25 2 + 2S. 



The second result follows from taking the square root, dividing by 



Qd\{w — d k ) 2 + 25 2 and bounding the remaining error terms. 
Remarks: 



□ 
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1. Note that in the case r\ = df., rf = 0, E = Dk and the integral in the second 
result above is divergent. We consider this case here. At Zq, we have 

w"(z ) = z r}/a,w"'(z ) = r}/a, 

and so there exists an ascent path for w" starting at z . By Pjj, this is also 
an ascent path for 2w 3 + 1. Following this path, even for a short distance 
(say of order e 1//14 ) allows us to reach a point z\, say, where we can use the 
above argument for small but nonzero w{z\) — dk- 

2. The above theorem extends easily to the case when rj ^ dk but rf = 0, 
because E — D k ^ 0. 

3. The case when rj = dk, but rj' ^ can be treated by an argument similar 
to the above ascent-path argument. 

4. The only remaining difficulty is that the path of integration may originate 
at or possibly intersect one of the simple zeroes of R(w). (Note that the 
definition of U does not exclude the simple zeroes p, a of R(w) from the 
path.) However, these problems are clearly overcome by ascent-path- type 
argument similarly to the way in which the zeroes of P(w) were avoided in 
Remark 4 following the proof of Theorem |2.1f 

These remarks cover all the remaining degenerate possibilities. 
3.3 Large Energy Limit 

Here we consider the only remaining case of possible limits of the generic 
solution, that is, we consider the case when \E\ (or equivalently r\ or rf, with 
if 7^ Tj j is large. 

Definition 3.6 Suppose e is an appropriate number. Parameters rj, rj 1 are 
called large if 

2E = T]' 2 - rf - 2r] 

satisfies 

\E\ > |lne| 4 . 

The corresponding solution w defined by initial values w(zq) = n w'(zq) = v( 
is called a large solution. 
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In this case, it is convenient to rescale variables as follows. Let 



w 
s 



E l ^q{s) 
E^z 



Then the first integral of Pjj becomes 

2 



9 s 2 = 9 4 + 2 + -^ I (g + J dssqq s /a). 



(3.32) 



Definition 3.7 Define the q-domain 

Q := {sE- 1 ^ eZ\\q\< \E\ l '\ \q A + 2| > \E\~ l l 4 } 
and s-domain 

§ := {s£~ 1/4 G Z | q(s) G Q}. 
e is called small enough if S is connected. 



Let s 



20 



E 1 / 4 , 



rjE 



-1/4 / 



r/i? 1 ' 2 . We will loosely call t, </, and 



= E~ 1 ^ 4 w(sE~ 1 ^ 4 ) large if 77, 77' and w(z) are large. 
Note that we will now restrict paths of integration to those with length 
upperbounded by 2nB in S. 

Theorem 3.5 Given an appropriate set {e,z ,B}, with |e/c|/2 < 1/2 where 
e is small enough, and large 1, d , and the corresponding solution q(s) satisfies 



dq 



so; 



< -|£|- 1/2 
2 

< 3nB\E\- 1/2 



for all s G S. Moreover, if Sq and sio are two successive points in S where 
w = 1, then for j = 1 or 2, 



Wj\ < 3ttB\E\- 1/2 . 



\{sw - s ) 

Proof: The small term in the first integral Eqn( |3.32D is given by 
T:=^ JI (q + sq 2 /(2a)- Jdsq 2 /{2c 



The results follow from the techniques of the proof of Theorem 2.1 



□ 
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A Complete Elliptic Integrals 

Here we prove some results for complete elliptic integrals as functions of E. 
Let C be one of Ck, k — 1, 2. 

Lemma A.l The elliptic integrals 

$ = j> dww 2 Vu> 4 + 2w + 2E, 



ip = <j> dw 



w 2 



* = -{ 



c Vw 4 + 2w + 2E' 

2 



dw 



Tc w 4 + 2w + 2E 3/2 ' 
considered as functions of E satisfy the following differential equations. 

d l = i, 
dE V 

*t = 1/ 
dE * 

d 2 i> ( 40E \ ~ 

IE* = [l28E*-27F (L33) 

Proof: The first two statements are obvious. To prove the last assertion we 
use the following definitions of other elliptic integrals. 



<j> dw Vu> 4 + 2w + 2E 



ic 

uo := & dw 



uj' := — <f> dw 



c Vw 4 + 2w + 2E 
1 



ic ^4 + 2^, + 2E 3/2 
:= jf dww Vw 4 + 2w + 2E 



IC 
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dw ■ 



c Vw 4 + 2w + 2E 
w 



dw 



<c w 4 + 2w + 2E 3/2 
That these are interrelated can be shown by integration by parts as follows: 

~ If, 3 2w 3 + 1 

ip = — f dww 



3 J c Jp(w 



2 r , 2 w A + 2w + 2E- (3/2)w-2E 
1 dww 2 



3 Jc .fpluj 



2 ~ r , w 3 4£, 
■-ty + f dw —j== + -z-V 

3 J c Jp w ) 3 



2~ l/,2^ 3 + l 1 4£ 

'o^+of dw /==== - + -S-V 
3 2Jc JP(w) 2 3 



Therefore, we get 

u = ^E1>- y^. (1.34) 
Similarly, starting with u and integrating by parts gives 

3£ = 30 + 4£cj, (1.35) 

and starting with gives 

20 =~%l) + 2E(j). (1.36) 
Finally, starting with 0, we get 

</>' = - (1.37) 
Differentiating Eqn( |1.34|) w.r.t. E yields 

J = j^V' - |V- (1-38) 
Doing the same with Eqn (|1.35|) gives 

- uj = 30' + 4Eu'. (1.39) 
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Now eliminating uo, u, uo' between Eqns( [Hl4* - 1.39 ) gives the desired differen- 
tial equation for ip(E). □ 

Note that these show that ujj and, therefore, the periods Uj are singular 
at the zeroes of 128-E 3 — 27. These are precisely the degenerate values D k 
found in Eqn( 2.12| ). 

Consider ipj, in the degenerate limit E — > D k (for some k). We drop 
the subscript j below, and consider momemtarily only the singular case. 
Frobenius expansion of the solutions ip of Eqn (|1.33|) show that there exists a 
singular solution with behaviour 

$ « co + a(E - D k ) \n{E - D k ) 

as E — > Dk, where Co is an arbitrary constant and c\ is a multiple of it. 
(For a derivation of similar results with more detail see Appendix B of |J.) 
Moreover, we have 

^ = -0'(£)~d \n(E-D k ) 
as E — >• D k . The period in this case is given by 

u = ±Ehl> - H 4, « l n (E _ (1.40) 

There is also a linearly dependent solution analytic at -Dfc with expansion 

ojj^ba + biiE-Dk) 
as — >• Dfc, where 6q is an arbitrary constant and 6i is a multiple of it. 
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